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In this lecture series I will give a fundamental insight into configuration space tech- 
niques which are of help to calculate a broad class of Feynman diagrams, the sunrise- 
type diagrams. Applications are shown along with basic concepts and techniques. 
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1 Introduction 



There has been a renewal of interest in calculating the loop integral of a simple topology, 
the so-called generalized sunrise (or sunset, water melon, banana, basketball) diagrams [H 
121 ini HI 1^ ■ The two- loop sunrise diagrams with different values of internal masses have 
been recently studied in some detail (see e.g. [HI 13 El El CHI IH] and references therein). In 
this lecture series I outline the configuration space technique IT! ^ IT U ITH l ITH l ITT t IT ^ ITil] . 
It can be used to verify the results obtained within other techniques [201 1211 1221 12H| and 
to investigate some general questions of Feynman diagram calculation [H [THl 1211 1211 1211 
123 I2H1 12ni IHOI- The diagram with the sunrise topology form a subset of the general set 
of diagram with a given number of loops [3T]. They appear in various specific physics 
applications: 

• sum rules for baryonic currents [S21ISSllSlllSSlESllS3ISHlEn] 

• gluonic correlators jlHl El 1121 

• multiloop calculations |131 C] 

• decays with multi-particle phase-space jH] 

• lattice QCD calculations |1S] 

• mixing of neutral mesons jlEj 

• Chiral perturbation theory (ChPT) and effective theories for Goldstone modes in 
higher orders of momentum expansion [Tfl E] 

• exotics |i8lll9] 

• effective potentials for symmetry breaking [HHl E] 

• finite temperature calculations jSH EH (SHI IHl IS21 ISS 

• applications in nuclear physics 

• applications in solid state physics [SH| 
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2 Concepts of configuration space techniques 



I will start this lecture series with explaining the basic elements that occur in calculating 
sunrise-type diagrams. The genuine sunrise diagram with three internal lines is shown on 
the left hand side in Fig. ^ It is the leading order perturbative correction to the lowest 
order propagator in 0^-theory. The corresponding leading order perturbative correction 
in 0'^-theory is a one-loop diagram and can be considered as a degenerate case of the 
prior example. A straightforward generalization of this topology is a correction to the free 
propagator in ^'^"'"^-theory that contains {n — 1) loops and n internal lines (see the right 
hand side of Fig. [T)) . 

mj 

// m2 ' \\ 
\\m„_j ■ / / 

Figure 1: Genuine sunrise and general topology of the class of sunrise- type diagrams 




In configuration space, the "calculation" of the correlator two-point function is very 
easy. It is just given by the product of the propagators, 

nix)=l[Dix,mi). (1) 

i 

The propagators themselves (for the scalar case and massive particles) are solutions of the 
Klein-Gordon equation and can be written in D-dimensional Euclidean space-time as 

^, , r d^p e*^^^^ (mx)^Kx(mx) im/xY , 

D(x,m) = / - — ^— = ^ , \, ' = ''. Kxinix) 2 

where I used D = 2\ + 2 for convenience, x = y/xj^, and Kx{z) is the modified Bessel 
function of the second kind (see e.g. ^\)- Frankly speaking, the propagator D{x,m) 
describes the propagation of a particle of mass m from the space-time point to x (a 
situation which of course is invariant with respect to translation). For m = we obtain 
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the massless propagator 

f dPp e*'^''^'' r(A) 

simply by calculating the limit, taking into account the series expansion of the Bessel 
function Kx{mx) which will be shown later on. Modified internal lines can also occur, as 
a result of a reduction procedure of more complicated diagrams or as diagrams with mass 
insertion on a line. This insertions are indicated by one or more dots placed on the line. 
Again, the representation of these cases is easily established, 

(note that D^'^\x,m) is related to the ji-ih derivative with respect to the parameter w?). 
It is obvious that, up to different power in x and normalization, a modification of the 
internal line just changes the index of the Bessel function. The modified propagator has 
the same functional form as the usual propagator. This matter of fact unifies and simplifies 
the considerations. 

The generalization to more complicated integrands with additional tensor structure 
pMi . . .p^^k jg straightforward as well. Because of 

d^p p^e^^^"^^ . d r d^p e^^^"^^ .Xi, d ( (mx)^ K x{mx)\ . , 



(271)-^ p2 + 777,2 Qx'' J (27r)^p2 + 777^ X dx \ {27r)^+^x'^^ 

this merely leads again to different orders of the Bessel function, together with a tensor 
structure x'^^ ■ ■ ■ x^''. However, these tensor structure does not matter in taking simply the 

product of the propagators. We can conclude that in order to obtain the expression for 
the correlator function in configuration space there is no integration at all. 

2.1 The correlator function in momentum space 

An integration occurs, if we want to calculate the correlator function in momentum space, 
being the Fourier transform of n(x), 

n(p) = j Ii{x)e'P'^'"'d^x. (6) 



For a general tensor structure, the integral is calculated by expanding the plane wave 
function exp(ip^x^) in a series of Gegenbauer polynomials C'^{p^x^ /px). The Gegenbauer 
polynomials are a generalizations of associated Legendre polynomials to the D-dimensional 
space-time. They are determined by the characteristic polynomial 



{e-2tz + l)-^ = Y.t'Ct{z) (7) 

1=0 

and are given by Cq{z) = 1, Ci{z) = 2Xz and the recursion formula 

(/ + 1)CU^) = 2{l + \)zC^{z) -{l + 2\- l)Cl,{z) = 0. (8) 
The Gegenbauer polynomials satisfy the orthogonality relations (xj = 

on the D-dimensional unit sphere with rotationally invariant measure d^X2 where 

Expressed in terms of Gegenbauer polynomials, the plane wave reads 

exp(^j9^a:^) = r(A) fP^y'f^^^{X + l)J^^,{px)C^{p^x>^/px) (11) 

^ ^ ^ 1=0 

where Jx{z) is the Bessel function of the first kind and p = y/PfiP^- This formula allows one 
to single out an irreducible tensorial structure from the angular integration in the Fourier 
integral. Integration techniques involving Gegenbauer polynomials for the computation of 
massless diagrams are described in Detail in jHO] where many useful relations can be found 
(see also jS3)- '^^^^ absence of explicit tensorial structure, however, the orthogonality 
relations for Gegenbauer polynomials can be used to explicitely integrate over the unit 
sphere, 

J e'P-^'d''x = 2n'+' (^^yj,{px). (12) 

In the following I will concentrate on this scalar case. The Fourier transform of the sunrise- 
type diagram, therefore, is given by the one- dimensional integral 

n(p) = 27r^+i jy (^^^ ^ Jx{px)D{x, mi) ■ ■ • D{x, mn)x^^^^dx. (13) 
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2.2 The renormalization 

While the configuration space expression for the sunrise-type diagram contains no integra- 
tion, it does not always represent an integrable function for a generic value of space-time 
dimension D. n(x) can have non-integrable singularities at small x for a sufficiently large 
number of propagators when D > 2 [Q2\. The reason is that each propagator by itself 
can be considered as a distribution. However, distributions do not form an algebra and 
multiplication is not well-defined. In this respect the correlator function n(x) is not com- 
pletely defined as a proper distribution. In attempting to integrate such a function over 
the whole x-space we encounter infinities which are usual UV divergencies. Therefore, 
the computation of its Fourier transform requires regularization (for instance, dimensional 
regularization) and subtraction. 

UV divergencies given as poles in e for the dimensional regularization are related to 
short distances x. Therefore, one should expand massive propagators at small x to obtain 
the counterterms. Still a convergence at large x should be retained which provides IR 
regularization (see e.g. 63 ). This technical constraint can be met for instance by keeping 
one (or two) massive propagators unexpanded. The possibility is given because the integrals 
of one or two Bessel functions are known IHH I(i5j . 

2 2 

However, one can do even simpler by using a damping function like e~'^ ^ to suppress 
contributions for large values of x, corresponding to IR singularities. This is done by intro- 

2 2 2 2 

ducing the factor 1 = e'^ ^ ^ into the integrand. In this case all massive propagators 

2 2 

along with the factor ^ can be expanded in small values of x. The expansion of the 
propagators is straightforward and results from the definitions 




(^/2) 



1 



{z/2f 



+ 



) 



(14) 



A;!r(A + A; + l) 



r(l + A) 



1 + A 



and 




(15) 



with 




(16) 
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Therefore, one has 

1 r(i-A)/^^2A- 



2J ^ ' 2 



1 + 



1 - A V2y r(i + A) V2 



+ 0(^^z'+'^). (17) 



The final integration can be done formally by using the identity 



" ^r'-ig-'^'^'fix - ^/x-T(r/2). (18) 

As an example we consider the gcnuin sunrise diagram with three different masses mi, m^-, 
and m3, given in configuration space by 

Yl{x) = D{x^m\)D{x^m'2)D{x^m-i). (19) 

The counterterms for the Fourier transform 

n(p) = 27r^+^ Jx{px)D{x, mi)D{x, m2)D{x, m3)x^^+Mx (20) 



,,2^2 ,,2^2 



are obtained by introducing e'^ ^ e ^ ^ and expanding the occuring Bessel functions as 

2 2 

well as ^ . The singular parts we obtain for dimensional regularization with D = 4 — 2e, 
i.e. A = 1 — £, are given by 

-4£ f ^2 I -Kv,2 I ^2 1 / 3 



. . _ At [ + mj + mj 1 / mf ln(mje^72//) \ 

iising W - 1 g^2£2 + £ l^Z. 128 1024 ) 



(21) 



where 7' = 7£;/2 — 3/4. Note, though, that the result is given in Euclidean domain. In 
this way one easily finds pole parts for any diagram with any mass arrangement. 

2.3 The momentum subtraction 

Renormalization by momentum subtraction is the oldest renormalization method. The 
idea is to subtract the integrand in some specific momentum point. For massive diagrams 
the point p = is IR-safe, and the receipe is realized for instance by expanding the function 

f)''^A(px) (22) 
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which is the kernel or weight function of the integral transformation) in a Taylor series 
around p = in terms of a p 
by writing 

and by keeping terms in the expansion on the right hand side. The number of 
necessary subtractions is determined by the divergence index of the diagram and can be 
found according to the standard rules of the /^-operation . For the example we started 
with, a momentum subtraction which is symmetric in the masses mi, m2, and ms leads to 
the finite part 

D{x, mi)D{x, m2)D{x, m3)a;^^+^e^'^''e"^''^'(ia; 

(24) 

where N is the order of expansion which has to be chosen for integrability. This is of 
course the same order which has been used earlier in order to extract the counterterms. It 
is obvious, though, that the expansion can cease at positive powers of x. Finally, e = 
can be taken to calculate the finite part. As a technical remark note that most numerical 
integration routines will run into problems if the upper limit is kept to be infinity. However, 
alread for values of the order x ~ 1 the integrand is negligible, so that the integration can 
terminate at this point. Also the region about the origin might cause trouble. In this case 
the integration interval can be subdivided close to the origin, and the whole integrand for 
the subinterval to zero can be expanded in a; as a whole. For specific values of p^, mi, m2, 
and ms we could reproduce results given in the hterature (e.g. results given in pT| IHE]). 

2.4 A second example 

The second example I want to show here starts with the same genuine sunrise diagram but 
is evaluated for a special setting of the masses. The reason is that there exists an analytical 
expression in the literature ^UJ- In this example the subtraction is done in several steps. 
The first step consists of a momentum subtraction for the weight function with A^ = 1. 
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olynomial series in p^. The subtraction of order A^ is achieved 



px\-^ ^ (-1)^ fpxY'^ 



nfin(p) = 27r 



A+1 



We obtain 



U{p) = n 

mom 

(p) + nrem(p) where 



(25) 



nmom(p) = 27r 



A+1 



2n 



A+l 



px 



px 



-A 



Jx{px) 



D{x, mi)D{x, m2)D{x, m3)x'^'^~^^dx 



J 1 



-A 



J\{px) 



1 p^x^ 1 

+ 



r(A + i) 4 r(A + 2) 



X D{x^mi)D(x,m2)D{x,ms)x'^^^^dx 



(26) 



nrem(p) 



A + p^B 

27^A+1 ^oo 



r(A + 



L I'OO 

/ D(x,mi)D(x,m2)D(x,ms)x'^^^^dx 

1 JO 



2^A+l ^oo 



/•OO 

/ x'^D{x,mi)D{x,m2)D{x,m3)x^^^^dx. (27) 

JO 



4r(A + 2) 

The analytical result in is given at the pseudothreshold p = mi + m2 — m^. For 
simplicity we choose mi = m2 = ^3/2 = m. In this case we obtain p = 0, nmom(O) = 
and the vanishing of the counterterm p^B. Only 

2^A+l ^00 



A = — / D{x,m)D{x,m)D{x,2m)x'^^~^^dx 

r(A + 1) JO 



(28) 



has to be determined. In order to separate singular and finite part as [271)^ A = S + F 
(where the total normalization of ^0] has been adopted) we use momentum subtraction 
for the last Bessel function K\{2mx) and obtain for the singular part 



S = ^^^^ J^ x'-''K,{mx)K,{mx 



r(A) 



1 + 



{mxf r(l-A), ,2A 



dx 



71 



4-2e 



(l_e)(l-2e) 
where we could make use of 



3 8 In 2 



— - + + 8(2 - 2 In 2 - In^ 2) 



+ 0(e) 



(29) 



/ x'^°'~^KJmx)KJ'mx)dx= — 7—— — -r(a + /i)r(a)r(a)r(a — /i). (30) 
JO m''°'T{2a) 



m2°r(2a) 

In comparing the part S with the analytical result in [lOj for p = 0, 



nref(O) = Tt'-'^ 



{l-e){l-2e) 



3 8 In 2 , 2 ' 
— - + 8 In^ 2 



+ 0(e) 



(31) 
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one obtains the same singular parts while the difference of the finite parts is expected to 
be found in the part F. Indeed, one obtains 



1) Jo 

r(A) 



r(A + 

X < {mx)^Kx{2mx) 



1 + 



A [mx ) X 

mxy Til-X)^_^,, 



1-A r(i + A) 



[mx] 



dx (32) 



which is finite and can therefore be calculated for = 4, 



F = 167rW/ —Ki{x)KAx)\xKi{2x)-- 1 + x^f-l + 27^ + 2 Inx) 
Jo X I 2 1- 



(33) 



It is shown numerically that F = 167r^m^(ln2 — 1), and because this quantity is finite, the 
normalization can be restored, we obtain 



F = vr 



4-2£ 



{l-e){l-2e) 



16(ln2 - 1) + 0(e). 



(34) 



2.5 A third example 

An example for the three-loop case is the sunrise-type diagram with two massive and two 
massless lines at vanishing external momentum. The analytical expression for the diagram 
in configuration space representation is given by 

2 / -rn/ , N \ 2 



n(0) = / {D{x,m)y{D{x,0)fd^x 



imx)^K\imx) 



r(A) 



d'^x. (35) 



(27r)^+ia;2^ / \4:Ti^+^x'^^ ^ 
While the angular integration in /^-dimensional space-time is trivial, 

27r^+i 



le'P^^'d^x = 2n'+'(^^y Jxipx) 



for p ^ 0, 



r(A + i) 

the problem of residual radial integration is solved by Eq. (jHUj) . The result reads 

n(o) = fa^Y'' 1 r(A)2r(i - A)r(i - 2A)2r(i - 3A) 



(36) 



(37) 



287r3A+3 r(A + l)r(2-4A) 

This result corresponds to the quantity Mi in which is the simplest basis element for 



the computation of massive three-loop diagrams in a general three-loop topology considered 
in [25 • We obtain agreement with 
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2.6 The spectral density 

From the physical point of view the interesting part of the analysis of sunrise-type diagrams 
is the construction of the spectral decomposition of the diagrams. For the two-point 
correlation function we determine the discontinuity across the physical cut in the complex 
plane of the squared momentum, = —m? ± iO. Taking for instance the n = 2 sunrise- 
type diagram with two different masses mi, m2 in arbitrary dimensions, in (Euclidean) 
momentum space given by 

r cf^k 1 

n(p) 



{2ti)D {k^ + ml){{k + pY + ml) 

r(l — A) /"V / X 2 / X 2 '2\^-^ 

[x[l — x)p — [1 — x)m^ + xm^j ax, (3J 



(4vr)A+i 

the physical cut has influence on the integral for x{l — x)p'^ -|- (1 — x)m'l + xm\ < 0. With 
p^ = —se^™ we can calculate the discontinuity to obtain 

Disc(x(l-x)p2 + (l-x)m? + xm^) = r(A)r(l - A) ^^^^ 

for xi < X < X2 where xi and X2 are the two zeros of the integrand which can be decomposed 
into s^'^{x2 — x)^^^{x — Xi)^^^. The spectral density is given by the discontinuity divided 
by 2Tii. Therefore, we finally obtain 

n(.\ - ^2A+i / {s - ml - ml)"^ - Amlml \ ^'^'^ 

^ ~ 4v/i(27r)2A+i 1^ As ) ^ ' 

where VL^ = 27i'^^'^ /T{d/2) is the volume of the unit sphere in rf- dimensional space-time. 
Eq. ()40j) is of need in the following which is the reason why I have given some details about 
the calculation. In the third lecture we will see that configuration space technique allows 
for some other efficient tool for determining the spectral density, namely by inverting the 
so-called i^-transform. In the next lecture I will continue with today's examples and give 
recurrence relations which allow to reduce the result to a finite set of master integrals, 
expressible in terms of transcendental numbers. 
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3 Recurrence relations and transcendental numbers 



I continue the previous lecture by introducing just three other examples for calculating an 

sunrise-type diagram with the help of configuration space techniques in order to compare 

them with results in the literature. A wide field of diagrams to compare with is a subclass 

of three-loop bubbles Bn [24 , namely 

f d^k d^l d^p 
BN{0,0,n3,n4,n5,n6) = / n/o-p/o n ^ 



m3^(7r^/2r(3-Z)/2)) 



^2?i3 m^"* m^'^s m'^"'^ 

with two propagators absent {rii = n2 = 0). Actually we choose the remaining indices 
77-3, -0.4, 77-5, and rig in a way that the results become finite, but we not only calculate the 
finite part but also the part porportional to e in order to compare with [21]. Written in 
configuration space, the particular subset of bubble diagrams I^at is given by 

POO 







It is obvious that we end up with integrals of products of four Bessel functions with non- 
integer indices and a non-integer power of x. The plan for today's lecture is as follows: 
After introducing the results for the three examples in terms of these integrals, I will 
show that at least the expansion up to first order in e is possible analytically. I will then 
introduce a reduction procedure in order to reduce the integrals to a basic set of two master 
integrals, 

U{r):= {K,{i)fCdi and L[{r) := {K,{x))' C He'^'^Cm dx (43) 
JO Jo 

and will add considerations on the general class of basic integrals Lnir) and iv^(r). 



3.1 The example 57v(0, 0, 2, 2, 2, 2) 

We start with the example which is represented by the central diagram in Fig. |21 Each of 
the lines is modified (indicated by the dots on the lines) which means that instead of the 
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Figure 2: The diagrams for Bn{0, 0, 2, 2, 2, 1), Bn{0, 0, 2, 2, 2, 2), and 5^(0, 0, 2, 3, 3, 4). 



propagators D{x,m) we have to use 

and obtain = ma;) 

We now can use the general formula [HS] 

to expand the Bessel function in a series with respect to its index. In case of K^ir{z), 
however, the first derivative vanishes and we obtain K_i;{z) = Kq{z) + 0{e'^). Therefore, 
in expanding 

(rn\\fro ^ = HMO)' e {l + {I + 2je - 2ln2 + 2lnOe + 0{e')) = 
{T{l + e)yT{2-e) ^ ' 

= 2{Ko{0)'e{l + {l + 2He^-^/2))e + 0{e')) (47) 

we obtain 

poo poo 

5^(0,0,2,2,2,2) = 2(1 + 6) / (iro(O) e'rfe + 4£ / (iro(0)dn(en/2)cie = 

Jo Jo 

= 2(1 + £)L4(3) +4£Li(3). (48) 
In comparing with the analytical result in [24j we checked numerically ^3] 

m = -^ + ^c(3). 
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3.2 The example Bn{0, 0, 2, 2, 2, 1) 

In the diagram on the left hand side of Fig. |21 one of the hnes is not modified. Therefore, 
we have to deal with one regular propagator factor 

DW(x,m) = Dix,m) = i^M^K.imx) = M^A^,_,(mx). (50) 

In this case we obtain 

B^iO, 0, 2, 2, 2, 1) = (r(i^^))3r(2-^) Jo K.^eiOe-'^'dt (51) 

Using Eq. pUj) . we obtain 

iri_,(0=i^i(0-^^o(0 + 0(5^) (52) 



and 

Ae{l + (l - 7 + ln(e^-e/2))e + 0(5^)). (53) 



(r(l + e))3r(2-5) ^ V +1 ^ 



The result reads 

E(0,0,2,2,2,l) = 4(1 + 5)/ {KoiOf K^{OedC - (Ko{0)CdC + 

Jo Jo 

+8e {M0fK,{0eHe'^^/m- (54) 
Jo 

This result is not yet written in terms of L^lr) and L''^{r) and will be treated after having 
introduced the reduction procedure. 

3.3 The example Bn{0, 0, 2, 3, 3, 4) 

In order to demonstrate the power of the configuration space technique also in this cal- 
culation, we finally choose the example found on the right hand side of Fig. |21 where the 
modified propagators 

D^'\x,m) = ^^^-i-e(mx), D^'\x,m) = ^^K_,_,imx) (55) 

are used, together with the expansions 



e_ ... _ ... _ . . 2e_ _ 2e 

Ko[x) + U{ 

(56) 
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We obtain 

5^(0,0,2,3,3,4) = s{r{l + e)m2 - e) i ^-(0 (^-i-.(O)' ^-2-.(0e'+'^c?e = 

1 _|_ (T /"CXD cr /"CXD 

+ 9670 ^^^^^^^^^^^^^'^'^^+4/0 (i^i(Ofe^rfe + 

+^ /o ^o(^) (^i(^)) ^2(0e'ln(e^-e/2)c^e- (57) 

3.4 The reduction procedure 

Especially in the last expression there are a lot of different integrals found which differ 
from the basis i^4(r) and L\{r). However, the integrands can be reduced to integrand in 
terms of -f^o(0 ^^'^ -^i(0 ^^h^ by using the relation 

71 — 1 

KniO = 2^X„_i(a;) + K^-2{x). (58) 

After the first step, namely the expansion of Bessel functions for non-integer indices, the 
above relation establishes the second step in our reduction procedure. Finally, we use 

^KoiO = -K,{0 and 

^^1(0 = -liMO + MO) = -MO-^^MO (59) 

to perform the third and last step, for instance 

Jo Jo 

roo poo 

= -3/ {Ko{OfK,{OCd^ + r {Ko{0)'r'd^ (60) 
Jo Jo 

and therefore 

L'i\r) = f (MOfMOCdC = \ f {KoiOfr'd^ = \U{r - 1). (61) 



16 



The complete set of reduction formulas is given by 





L^\r) = 


{r-3)Lf\r- 


l)-3Lf (r), 






L?\r) = 




r-l)-2Lf{r)] 






L?\r) = 




r - 1) - U{r)) , 






L^\r) = 


I 








Lf\r) = 


{r-3)Lf\r- 


+ (r-1) 


-3Lf (r), 




Lf\r) = 


^((r-2)Lf 


(r-l) + Lf (r- 


1) - 2Lf ) 




Lf\r) = 




{r-l) + nl\r- 


1) - L\{rX 




Lf\r) = 


I {rL\{r - 1) - 


f ^4(r-l)), 




in general 










LtHr) = 
LT\r) = 


' ( 


(r-m + l)4" 
(r-m + l)L^[(™ 


'^\r — 1) — (m — 
-i)(r-l) + 4'"- 


i)(r- 1) - 


n — m + 1 ^ 

' I 


n — m + I ^ 



(62) 



(63) 



(64) 



and is coded in MATHEMATICA in order to automatically reduce to the master inte- 
grals pni- The results for the examples after executing the second and third step of the 
recursion reads 



5^(0,0,2,2,2,2) = 2{l + e)L4{3)+4eL\{3) + 0{e^ 

5^(0,0,2,2,2,1) = 2Li{l)+4eL[{l) + 0{e^), 
5iv(0,0,2,3,3,4) -- 



^^.(3)-^L,(5)-^L,(7) + 0(.). 



Comparing with the results obtained by using the RECURSOR package 

5^^(0,0,2,2,2,2) = + 1({3) + 0{e), 

8 16 

5;v(0,0,2,2,2,l) = l({3) + 0{e), 

5;v(0,0,2,3,3,4) = ^ + 0{e) 

576 



(65) 



(66) 
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(the order e is omitted here) we can adjust the master integrals 

L4(3) = + T^C(3), ^4(1) = ^C(3), 16L4(3) - 2^4(5) - = 1. (67) 

lo 32 8 

Especially the last relation is interesting because it means that 

roo 1 

/ K,{i) {K,{i)f K,m' di = -. (68) 
JO 3 

This surprizing identity has been checked numerically as well. Note that only Li{r) for odd 
values of r are needed for the leading order contribution. However, before concentrating on 
the master integrals themselves, let me add some remarks on the efficiency of the reduction. 

3.5 The efficiency of the reduction 

Taking to be the maximal total power of denominator factors, the calculation of a three- 
loop integral like in the previous examples would create and deserve a field of integrals 
which are recursively defined by each other using recurrence relations [21]. In contrast to 
this, the reduction procedure in our case needs steps to reduce the result to a basis of 
2[A^/2] — 5 basic elements L^ir) and L\{r). Therefore, the "costs" of our method are of 
the order A^^ which considerably reduces the time consumption in a computer evaluation. 

3.6 Values for the master integrals 

By comparing the results for different sun-rise type diagrams obtained by using our reduc- 
tion method [T3] and the procedure RECURSOR in |21] we can find values for the master 
integrals. The first results read 

7 7 3 49 27 

^4(1) = ^C(3), ^4(3) = ^C(3) - ^, L4(5) = ^C(3) - • • • (69) 

In general we found that the (odd-valued) master integrals are given by expressions of 
the kind L4(2r + 1) = Aj.C{'i) — B^. Of course we might ask whether there is a general 
way to calculate these master integrals. At least we can execute the calculation for the 
master L/^{1). This master integral is (by a fresh view, i.e. not considering where the 

18 



master integral comes from) given by the three-loop sunrise-type bubble in two space-time 
dimensions (i.e. A = 0), 

n(0) iself can be written as the two-dimensional integral over the product of two two-line 
correlators of the kind 

2ny/j?\/p^ + Arri^ \ y/p^ + 4m^ — ^/p^ ) 47rm2sinhr 
(the last identity for substituting p = 2m sinh(r/2)), so that 

^ ,2.).W/^ (n.„))^ -\r^- 5C(3) m 

where jHlI 

r"-Mr 2° - 1 , , , 

sinhr 2°-i ^ ^ ^ ^ 

has been used. The expression for general (odd) values for r reads 

L4(2r + 1) = 27im'^ J Il{p) (^-m^npY Il2ip)d^p (74) 

where Dp = d"^ /dp^dp^ is the two-dimensional d'Alembert operator in (Euclidean) momen- 
tum space. It is still an open question whether this integral can be calculated to obtain the 
simple expression containing rationals and C(3)- However, a closer look at the numerical 
value shows that, apart from the chance to obtain an analytical result, the attempt might 
not be worth the effort. Indeed, looking at the master integrals for higher and higher values 
of r we find that the contributions y4rC(3) and Br increase while the difference leads to a 
significant cancellation, for /4(11) for instance up to three significant figures. Instead of 
trying to calculate the analytical expression, we therefore look for an asymptotic expansion 
of the integral at large values of r, 

i*(2''+ 1) = (1-73^75 + 0(1/'-')) • (75) 

We can introduce a new parameter k which accounts for higher order terms in the asymp- 
totic expansion. The expression 

, 7r2r(2r) / 1 \ 

i^(2'-+i) = ^(i + ;T^) (76) 
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with K = 0.97 gives numerical results with an accuracy better than 1% for all r > 1, while 
for r > 3 the relative accuracy is better than 10^'^. A similar expression can be found for 
the logarithmic masters, 

L[{2r + 1) = ^ m2r) + - 3 ln2) (l - (77) 

where "^{z) = T'{z) /T{z) is the polygamma function or logarithmic dericative of the gamma 
function. For ni = 1.17 the accuracy for r > 3 is better than 1%, for r > 5 better than 
10~3. 



3.7 Back to the genuine sunset 

As it is already visible in Eq. ()64p. the recurrence relations just discussed for n = 4 massive 
internal lines are given for other numbers of internal lines as well with no costs. We could 
for instance return to the genuine sunset with n = 3. In this case the master integral -^^3(1) 
is given by a two-loop sunset-type bubble in two space-time dimensions, 

Ti\ + l)Jo [ (27r)^+ix2A ) '^'^ (27r)2m2yo ^ ^^^^ 

and, therefore (t = e""^), 

r N 2 / ^2{p) ,2 f"^ rn^rdr \ntdt , 

L,il;M/m)=m I ^^,d p = 4^. ,i,,.(,/2) + = " io 1 - 27t + ^''^ 

where the mass M of the third line is kept different from the other two masses m for a while. 

By differentiating with respect to M, any positive power of the propagator (and/or power 

in X in configuration space) can be obtained - therefore, all master integrals. Because of 

7 = 1 — M^/2m^, the denominator has two roots ti^2 = 7 ± y/'j'^ — 1 with tit2 = 1, and we 

obtain 

Ul;M/m) = -r '^'f n.(l/tO-n.(l/t2)^Li2(t.)-Li2(t.) 

where Li2(2;) is the dilogarithm function. The differentiation with respect to M is now 
straightforward and can be performed with a symbolic manipulation program. The pseudo- 
threshold case M = 2m is a degenerate case, we obtain 
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while for M = m we obtain 7 = 1/2, ti^2 = exp(±i7r/3), and therefore 



Ul) = LM- 11 



1 \ntdt 



/o l-t + t2 73 
a result in terms of Clausen's dilogarithms. 



^^ImLi2(e 



f82) 



3.8 Generalization to the spectacle topology 





Figure 3: "spectacle+propagator" 
representation, also called "specta- 
cle" topology diagram, in two differ- 
ent forms. The configuration space 
points 0, X and y are indicated. 



To finish today's lecture, we can look a bit beyond 
the rim of the plate of sunrise-type diagrams. Ac- 
tually, we can generalize the configuration space 
technique to the so-called spectacle topology, as 
shown in Fig. 3 in two different representations. 
The configuration space expression of a spectacle 
topology diagram (again with a different mass M 
on the frame) written in a form suitable for the 
actual purpose is given by 



S{M) = J D{x-y,M)D{x,mfD{y,mfd^xd^y. 



(83) 



The key relation for a significant simplification of the configuration space integral with 
spectacle topology is given by a relation obtained from the addition theorem of Bessel 
functions jHS]- I will not go into detail here, only mention the result ^H] 



^A(|r-p|) _ h{p)K^{r) 
(27r)^+i|r-p|^ ' 



r > p 



where I\(z) is the modified Bessel function of the first kind. This result allows one to write 



S{M) 



M^^ Kx{M\x-y\) 



D{x, mfD{y, mfd^x d^y 



(27r)^+i {M\x-y\)^ 
' 27r^+i \^ /-oc 



{p{y,m)f x^^^^dy x 



{MyY {MxY 
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As an illustration we look at the two-dimensional space-time where 

1 /-oo 2 f°° 2 

S{M) = j—:^ x{Ko{mx)) dx y {Ko{my)) dy x 

X {KoiMx)Io{My)e{x - y) + Ko{My)h{Mx)9{y - x)) , (86) 

comparing this result with the momentum space representation 

^ ^ J (27r)2 p2 + m^io (1 -t2)(l -27t + t2) ^ ' 

where the same chai of substitutions have been used. We obtain 

S(M,^/M^, /W . ™Z|^(f±l)M (88, 

Ti — 12 TTL [t — ij 

where Lis (2;) is the trilogarithm function. While for M = m we obtain a result including 
the Clausen trilogarithm Cl3(27r/3), for the case M = 2m the integral for S{M = 2m) 
simplifies as in the case of the sunrise diagram, one obtains 

'5(2m) = £jgc(3)-ln2). (89) 

In all these cases the sixth order roots of unity ±1, exp(±7r/3), and exp(±27r/3) play an 
important role [HZ]. In P^l^ we have followed the lines given in and have constructed 
a shuffle algebra for the occuring integrals which finally lead to the transcendental numbers 
which normally occur in these calculations. At this point I will not go into detail. Instead I 
close the lecture by mentioning that the numerical comparison between the result obtained 
in the momentum space approach with the result within the configuration space technique 
show conicidence. Tomorrow we will deal with the consideration of sunrise-type diagrams 
close to threshold. 
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4 Expansions close to threshold 



From the physical point of view the interesting part of our analysis of sunrise-type diagrams 
is the construction of the spectral decomposition of the diagram. Using the configuration 
space technique, the spectral density can be constructed directly from the correlator func- 
tion in configuration space. The technique introduced in is based on an integral 
transform. Starting with the dispersion relation in (Euclidean) momentum space, 

~ = r 4^ = r p{s)b{q, V^)ds, (90) 
Jo + s Jo 

the corresponding relation in configuration space reads 

roo 

U(x) = / p{s)D{x,^/s)ds. (91) 
Jo 

This representation was used for sum rule applications in [^11113 where the spectral density 
for the two-loop sunrise diagram was found in two-dimensional space-time [70]. With the 
explicit form of the propagator in configuration space given by 

D{x,V~s)=^-0^Kx{V~sx), (92) 

the representation turns into a particular example of the Hankel transform, namely the 
i^'-transform [TH [721 

9{y) = / f{x)K^{xy)y/xydx. (93) 
Jo 

The inverse of this transform is known to be given by 

]^ rc+ioo 

f{x) = — g{y)h{xy)y/xydy (94) 

ITT J c—ioo 

where the vertical contour in the complex plane is placed to the right of the right-most 
singularity of the function g{y) [72j. Translating this inverse i^-transform to our situation 
which relates p(s) and n(a;), we obtain 

, pc+ioo 

p{s) = -(27r)\ / U{x)x^+^h{^/sx)dx. (95) 

J c—ioo 

The inverse transform given in Eq. (j95p solves the problem of determining the spectral 
density of sunrise-type diagrams by reducing it to the computation of a one-dimensional 
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integral for the general class of sunrise- type diagrams with any number of internal lines and 
different masses. Because the contour can bypass the area of small values of x, the integral 
is finite. Using this expression, the spectral density is known numerically as a function 
of s. Eq. (j^J3j) will be used in the following to compare the exact result with different types 
of expansions found in the literature, in order to test their convergence behaviour. 

4.1 Considerations in Minkowskian domain 

The threshold region of a sunrise- type diagram is determined by the condition g^ + M^ = 0, 
where q is the Euclidean momentum and M = is the threshold value for the spectral 

density. We introduce the Minkowskian momentum p defined by p"^ = —q^ which is an 
analytic continuation to the physical cut. Operationally, this analytic continuation can 
be performed by replacing q —>■ ip. To analyze the region near the threshold we use the 
parameter A = M —p which takes complex values, while in phenomenological applications 
the parameter E = —A=p — Mis more convenient. The spectral density will be written 
as a function of E in the following, p{E) = p{{M + E^). The analytic continuation of the 
Fourier transform Il{q) to the Minkowskian domain has the form 

Il^{p) = 27r^+i 1^°^ (^^) ^ Jx{ipx)U{x)x^^+^dx (96) 

(as we work in Minkowskian domain only, the index M will be dropped in the following). 

4.2 Large x behaviour of the weight 

For the threshold expansion we have to analyze the large x behaviour of the integrand. It 
is this region that saturates the integral in the limit p ^ M or, equivalently, E ^ 0. It 
is convenient to perform the analysis in a basis where the integrand has a simple large x 
behaviour. The most important part of the integrand is the Bessel function J\{ipx) which, 
however, contains both rising and falling branches at large x. It resembles the situation 
with the elementary trigonometric function cos (2;) to which the Bessel function Jx{z) is 
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rather close in a certain sense, as we will see. cos(2;) is a linear combination of exponentials, 

cos(^) = ^ (e'' + e-'') . (97) 

In the same manner, the Bessel function Jx{z) can be written as a sum of two Hankel 
functions H^{z) — J\{z) ± iY\{z), 

Uz)^^{hUz) + H^(z)), (98) 

and for pure imaginary argument the Hankel functions show simple asymptotic behaviour 

H^{iz) - z-^^^e^\ (99) 

Accordingly, we split up n(p) into n(p) = Il'^{p) + n~(p) with 

n±(p) = TT^+i (^^) ^ H^{ipx)U{x)x^^+'dx. (100) 

The two parts n^(p) of the polarization function n(p) have completely different behaviour 
near threshold and are analyzed independently in the following. 

4.3 The polarization function 11^ (p) 

The behaviour of Il~^{p) for large x is given by the asymptotic form of the functions which 
I simply write up to the leading terms as 



H+(ipx) = W^e-f^(l + 0(x-^)), Kimx) = .[^—e-"'^(l + 0(x-^)). (101) 
V nrpx V zmx 

We conclude that the large x range of the integral (above a reasonably large cutoff param- 
eter A) has the general form 

nX(M-A)~/ x-^e-^^^-^^'^dx, a = (n-l)(A+l/2). (102) 

J A 

The right hand side is an analytic function in A in the vicinity of A = 0. It exhibits no cut 
or other singularities near the threshold and therefore does not contribute to the spectral 
density. I could therefore go ahead to the second contribution n~(p). However, I want 
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to do a bit more in order to fill n+(p) with "life" (i.e., meaning). In using the relation 




(103) 



Figure 4: Representation of the reg- 
ular part n+(— p^) as vacuum bubble 
with added line. The cross denotes 
an arbitrary number of derivatives of 
the specified fine. 



between Bessel functions fo different kind I can 
replace the Hankel function H^(ipx) by the Bessel 
function Kx{px). But since the propagator of a 
massive particle (massive line in the diagram) is 
given by the Bessel function K\{mx) up to a power 
X, the weight function behaves like a propagator 
of an additional line with "mass" p. The explicit 
expression is given by (cf. Fig. 4) 



n+(p) 



-2Tii 



i2A+l 



(p2); 



n+(x)x2^+Ma;, n+(a;) = Tl{x)D{x,p). 



(104) 



4.4 The polarization function 11 (p) 



In contrast to the previous case, the integrand of the part 11 (p) contains H [ipx) which 
behaves like a rising exponential function at large 



H-{ipx) ~ x-^'^e^"". 
Therefore, the integral has the large x behaviour 

nx(M- A) ~ r x-^e-^ 

J A 



(105) 



dx. 



(106) 



The exponential factor e~^^ which estabhshes an effective inverse proportionality between 
A and x is the reason for the importance of the large x behaviour as being essential for 
the near threshold expansion of the spectral density. However, for A < the integral 
diverges at the upper limit, leading to non-analyticity at the threshold A = 0. In the 
complex A plane with a cut along the negative real axis the integral is analytic. This cut 
corresponds to the physical positive energy cut. The discontinuity across the cut gives rise 
to the non- vanishing spectral density of the contribution II~ (p) . 
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In order to obtain an expansion for the spectral density near the threshold in an an- 
alj^ical form we make use of the asymptotic series expansion for the function 11 (x) that 
crucially simplifies the integrands but still preserves the singular structure of the integral 
in terms of the variable A. Because of the already mentioned asymptotic expansion 



K^{mx) - V ^j^e-"*^ (l + 0{x-^)) , (107) 

the asymptotic expansion of the function n(a;) consists of an exponential factor e^^^ and 
an inverse power series in x up to an order N closely related to A'". It is this asymptotic 
expansion that determines the singularity structure of the integral. We write the whole 
integral in the form of the sum of two terms again, 

n-(p) = 7r^+i|(^!|^) '//,-(ipx)(n(x)-n-(x))x2^+^+^^dx + 

+7r^+' / (^) ^ H^{ipx)Ii''N{x)x^^+^+^'dx = Ii''\p) + II"^(p). (108) 

The integrand of the difference term 11°'* (p) behaves as 1/ x^ at large x while the integrand 
of the second term accumulates all lower powers of the large x expansion. We will again 
consider the two parts in turn. However, a comment on the regularization is in order here. 

4.5 Comment on the regularization used 

Being interested only in the large x behaviour, we introduced a cutoff A which necessary 
at least at this point because the asymptotic expansions are not defined for small values 
of x and lead to divergences. However, from the practical point of view the calculation of 
the regularized integrals with an explicit cutoff is inconvenient. Instead, the standar way 
to cope with such a situation is to use dimensional regularization. Note that dimensional 
regularization does not necessarily regularize all divergences in this case (in contrast to the 
standard case of ultraviolett divergences) but nevertheless is sufficient for our purposes. 
Therefore, we use a parameter e independent of A to regularize the divergences at small x. 
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4.6 The difference part n^^(p) 

In n*(p) the subtracted asymptotic series to order N cancels the inverse power behaviour 
of the integrand to this degree N. Therefore, the integrand decreases sufficently fast for 
large values of x and the integral even converges for A = 0. The part n*(p) is regular and 

gives no contribution to the spectral density up to the order A^. 

4.7 The asymptotic part n"^(p) 

The integral 11"* (p) is still rather complicated to compute. But we can go a step further 
in its analytical evaluation. Indeed, since the singular behaviour of 11"* (p) is determined 
by the behaviour at large x, we can replace the first factor, i.e. the Hankel function, in the 
large x region by its asymptotic expansion up to some order N, 



2 




with 



(-l)"(A,n) J-inX,N) 



(109) 



to obtain the double asymptotic representation 

n'^»*(p) = 7r^+^ (^^) ^ H^;j^{ipx)U''j^ix)x''^+^+^'dx. (Ill) 

Both asymptotic expansions are straightforward and can be obtained from standard hand- 
books on Bessel functions. We therefore arrive at our final result: the integration necessary 
for evaluating the near threshold expansion of the sunrise-type diagrams reduces to inte- 
grals of Euler's Gamma function type, i.e. integrals containing exponentials and powers 
in the integrand. Indeed, the result of the asymptotic expansion of the integrand is an 
exponential function e~^^ times a power series in 1/x, namely 

N-l A . 
a+2e^-Ax \ " 



X ' e 



T. -f (112) 



where a = (n — 1)(A + 1/2) has already been used earlier and the coefficients Aj are 
simple functions of the momentum p and the masses m^. This expression can be integrated 
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analytically using 



roo 

/ x-''+^'e-^''dx = T{l-a + 2e)A^-^-^'. (113) 

^0 



The result is 

N-l 

Yldas^^ _ A) = ^ AjV{l -a-j + 2£)A'^+^-^-2^ (114) 

j=0 

This expression is our final representation for the part of the polarization function of a 
sunrise-type diagram necessary for the calculation of the spectral density near the produc- 
tion threshold. Starting from this main result in Eq. (jll4p . I discuss the general structure 
in some detail. In the case where a takes integer values, these coefficients result in 1/e- 
divergences for small values of e. The powers of A in Eq. pi4p have to be expanded to 
first order in e and give 

— A^^ = — + lnA + 0(£). (115) 

Because of 

Disc ln(A) = ln(-E - iO) - \n{-E + iO) = -27iie{E) (116) 

Il^°-^[M — A) in Eq. (jll4j) contributes to the spectral density. For half- integer values of a 
the power of A itself has a cut even for e = 0. The discontinuity is then given by 

DiscVA = -2iVE9{E). (117) 

Our method to construct a threshold expansion thus simply reduces to the analytical 
calculation of the part 11'^"* (p) which can be done for arbitrary dimension and an arbitrary 
number of lines with different masses. In the next subsections I work out some specific 
examples which demonstrate both the simplicity and efficiency of our method. 

4.8 First example: equal mass (n = 3) sumise diagram 

For the genuine sunrise diagram with equal masses (the last restriction for reasons of 
representation only), the double asymptotic expansion reads 



-1 ™3/2p(p-3m)a; 

-M^) H-,ipx)U-ix)x^-^^ = ( J^3^3/. X 



1+ ^ 3^ 9 27 15 _^Q^^-3A Q^gN 

8mx 8px 128m^a;^ GAmpx^ 128p^x^ J 
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exact function 
2nd order expansion 
3rd order expansion 
12th order expansion 



0.2 



in units 1 0" 



0.1 



3 m, D = 4 



0.2 0.4 0.6 0.8 1 

7] = E/M 

Figure 5: Various results for the spectral density for n = 3 equal masses in D = 4 space- 
time dimensions in dependence on the threshold parameter E/M. Shown are the exact 
solution obtained by using Eq. (solid curve) and threshold expansions for different 
orders taken from Eq. pi9|) (dashed to dotted curves). 

The spectral density is obtained by performing the term-by-term integration of the series 
and by evaluating the discontinuity across the cut along the positive energy axis E > 0. 
The result reads 

f. 1 7 2 3 3 39 4 57 5 
"'^^ = SSVll' - 2" + 16" - 8" + I28" - 256" <"^> 

1024' 256' 32768' 65536' 524288' ^' 

where the notation rj = E/M, M = 3m is used ^]. The simplicity of this derivation is 
striking and easily reproduces the expansion coefficients in [10 . For different orders the 
result is compared with the exact result from inverse iiT-transform (Eq. ()95p) in Fig. El 
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Figure 6: Various results for the spectral density for n = 4 equal masses in D = 4 space- 
time dimensions in dependence on the threshold parameter E/M. Shown are the exact 
solution obtained by using Eq. (solid curve) and threshold expansions for different 
orders taken from Eq. p20|) (dashed to dotted curves). 

4.9 Second example: equal mass n = 4 sunrise-type diagram 

The sunrise-type diagram with four or more propagators cannot be easily done by using 
the momentum space technique because it requires the multiloop integration of entangled 
momenta. Within the configuration space technique the generalization to any number of 
lines (or loops) is immediate by no effort. For the spectral density of the equal mass n = 4 
sunrise-type diagram we obtain (cf. Fig. ^ 

EV^M'^/^ 1 81 o 2811 , 17581 , 



^ ... , . 81 2 2811 o 

piE) = ^{1 71 -\ ri^ ri-^ + 

268807r5y2l 4 ' 352 ' 18304 ' 292864 



V 



(120) 



1085791 , 597243189 . 4581732455 7 496039631453 „ 



19914752 ' 3027042304 ' 12108169216 ' 810146594816 
where t] = E/M and M = 4m is the threshold value. One sees the difference with the 
previous three-line case. In Eq. ()12()|1 the cut represents the square root branch while in 
the three-line case it was a logarithmic cut. This underlines my earlier considerations. 
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4.10 Strongly asymmetric mass arrangement 

As we have seen, the threshold expansion for equal (or close) masses breaks down for 
E ^ M = rrii. However, if the masses are not equal, the region of the breakdown of 
the expansion is determined by the mass with the smallest numerical value. Actually, one 
finds a breakdown at ~ 2mo where mo is a mass much smaller than the others. To our 
best knowledge, this issue has not been touched earlier. In the following I will describe a 
technique called resummation of the smallest mass contribution Instead of the double 
asymptotic expansion we use the partial asymptotic expansion 

U^^^ip) = ^'^'j^ (^) \-,;^^{tpx)IiZ,{x)x''^'^'^dx (121) 

where the asymptotic expansions are substituted for all propagators except for the one 
with the smallest mass mo, 

n-(x) = n-i(x)D(mo,x). (122) 

Because of the knowledge of 
Jo 

- 2 '^+2'^~5^j ^^^^^ 

with a = A — mo and (3 = rriQ, 11^"* (p) is expressible in terms of hypergeo metric func- 
tions [75j. The same is valid for the spectral density, because 

Disc / x''-^e''''K^{px)dx = 

2Txi Jo 

«V2-./5. T{3/2-fx)''[ 2 ' 2 ' 2 ^' j l^^^^ 

where a = E + mo. For the case of an n = 2 sunrise-type diagram with masses m and 
mo <^ m we compare the exact result in Eq. for the spectral density with the pure 
expansion near threshold 



-das 



-aasrr,. ^ V2^^^ i, , fl I 7\E 
^ ^ ' M^M^I^ I Vm mo M; 4 



M3/2 \ ^ \m mo Mj 4 
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Figure 7: Various solutions for the spectral density for two masses m and mo <^ m and 
D = A space-time dimensions. Shown are the exact solution which is obtained by using 
Eq. (solid curve), the pure threshold expansions using Eq. ()125|) (dotted curves), and 
the solutions for the resummation of the smallest mass contributions like in Eq. ()126|) 
(dashed curves), both expansions from the first up to the fourth order in the asymptotic 
expansion. For the pure threshold expansion the order is indicated explicitly. 

(M = m + mo, the second order asymptotic expansion should suffice to show the general 
features in a short and concise form) with the resummed result 



mE{E + 2mo) f / 1 3 



m: 



E{E + 2mo) (1 5 ml \ 

+ 8m(E + M) I 2' ^' 2' ^ " (E + mo)2 i ^^^^^ 



2 



E\E + 2moy / , lQm{E + M)\ ^ / 3 7 ^ m^ 
+ ^TTT- ^1 2-ri|i,-;-;i 



128m2(E + M)2\ 5(E + mo)2 / V'2'2' (E + mo) V J 

Having used e = for the finite spectral density, the first term in curly braces of Eq. 
is obviously equal to 1 in this limit because the first parameter of the hypergeometric 
function vanishes. However, I keep Eq. ()126p in its given form to show the structure of 
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Figure 10: The spectral density for the four-hne sunrise-type diagram with equal masses 
for D = e = 2.718 . . ., D = 3, and = vr = 3.14 . . . space-time dimensions, to demonstrate 
the practical convenience of the method. 

the contributions. The result of the comparison is shown in Fig. [7| The convergence of 
different orders of resummed expansions is remarkable, while the pure exansion diverges 
already at an early stage. 

4.11 A last diagram 

The last figure I want to show is a funny one. In order to demonstrate the universality 
of the configuration space method, in Fig. I show the spectral density for the n = 4 
sunrise-type diagram obtained via the inverse fT-transform in Eq. (j95j) for the space-time 
dimensions D = e = 2.718 . . D = 3, and D = tt = 3.14 . . .. 
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